Eccentricity fluctuations in an integrated hybrid approach: Influence on elliptic flow 
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The effects of initial state fluctuations on elliptic flow are investigated within a (3+1 )d Boltzmann 
+ hydrodynamics transport approach. The spatial eccentricity (erp and e part ) is calculated for 
initial conditions generated by a hadronic transport approach (UrQMD). Elliptic flow results as a 
function of impact parameter, beam energy and transverse momentum for two different equations 
of state and for averaged initial conditions or a full event-by-event setup are presented. These 
investigations allow the conclusion that in mid-central (b = 5 — 9 fm) heavy ion collisions the final 
elliptic flow is independent of the initial state fluctuations and the equation of state. Furthermore, 
it is demonstrated that most of the V2 is build up during the hydrodynamic stage of the evolution. 
Therefore, the use of averaged initial profiles does not contribute to the uncertainties of the extraction 
of transport properties of hot and dense QCD matter based on viscous hydrodynamic calculations. 
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I. INTRODUCTION 

After the observations at the Relativistic Heavy Ion 
Collider (RHIC) that the elliptic flow measurements are 
consistent with ideal fluid dynamics predictions the 
development of numerical solutions of viscous hydrody- 
namics equations is actively pursued The idea of 
those calculations is to quantify the deviations from local 
thermal equilibrium by comparisons to the available data 
for spectra and elliptic flow. 

Solving the differential equations of any kind of fluid 
dynamics prescription implies the knowledge of the 
boundary conditions, i.e. the initial conditions and the 
freeze-out criterion. Therefore, it is essential to investi- 
gate in a systematic way how different initial profiles and 
freeze-out implementations influence the final observable 
results. This article focuses on this question concerning 
the initial conditions. 

The realistic situation of the collision of two nuclei suf- 
fers from many different sources of initial state fluctua- 
tions @. The density profiles are not smooth, but there 
are peaks in the transverse and the longitudinal direction. 
There are impact parameter fluctuations within one spe- 
cific centrality class leading to multiplicity fluctuations 
and differences in the initial geometry 10]. Furthermore, 
the nuclei do not necessarily collide in the event-plane 
given by the laboratory system, but might also have a 
rotated reaction plane. All these effects arc averaged 
out, if assuming a smooth symmetric initial density pro- 
file as it is widely done by parametrising the initial con- 
ditions for hydrodynamic calculations within a Glauber 
or Color Glass Condensate (CGC) picture. To avoid 
potential misunderstandings, Monte-Carlo Glauber ap- 
proaches (and CGC approaches) that produce fluctuating 
initial conditions are available and widely used, e.g. by 
experimental collaborations. However, usually these fluc- 
tuating initial conditions are not used as event-by-event 
initial conditions for the majority of hydrodynamic sim- 
ulations. For alternative hydrodynamic approaches with 



fluctuating initial conditions see e.g. [Til H^. 

In this article an integrated Boltz- 
mann+hydrodynamics transport approach is applied 
to the simulation of heavy ion reactions in the energy 
regime from E\ a b = 2 — 160A GeV. The initial condi- 
tions are generated by the Ultra-relativistic Quantum 



Molecular Dynamics (UrQMD) approach [13|, |14j and 
the above mentioned event- by- event fluctuations are 
further propagated in the ideal hydrodynamic evolution 
employed for the hot and dense stage of the collision. 
On the other hand, calculations with averaged initial 
conditions are performed using the very same general 
setup. Comparing these calculations to the fluctuating 
setup the effect on the final observable elliptic flow is 
estimated. 



II. INITIAL ECCENTRICITY 

Let us start by looking at the initial state eccentricity 
as it is given by the UrQMD approach at the starting time 
^start = 2i?/ ^/■f 2 — 1 of the hydrodynamic evolution. For 
the specific values of these times for Au+Au/Pb+Pb col- 
lisions at different beam energies see Table HI This is the 
earliest possible transition time at which local equilib- 
rium might have been established after the two nuclei 
have passed through each other. 

The eccentricity quantifies the spatial anisotropy of 
the initial state which is transformed via pressure gra- 
dients into a momentum space anisotropy in the trans- 
verse plane that can be quantified by the elliptic flow 
coefficient. The standard definition for the eccentricity 
is the reaction plane eccentricity 



(1) 



with a 2 — (x 2 ) - (x) 2 and a 2 = (y 2 ) ~ (y) 2 . Especially 
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12.417 
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TABLE I: Starting times of the hydrodynamic evolution for 
Au+Au/Pb+Pb collisions at different beam energies. The 
eccentricities displayed in Fig. [T] are also calculated at these 
times. 



for smaller colliding systems the so called participant ec- 
centricity is popular 
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where an additional correlation term a xy — (xy) — 
(x)(y) is introduced. The participant eccentricity incor- 
porates the fact that the participants themselves might 
be rotated with respect to the reaction plane. The aver- 
ages (•) indicate averages over all particles in one event 
and the mean eccentricity is obtained by averages over 
many UrQMD events. The standard deviation (shown as 
error bars in FigJTJ) with respect to the event average is 
5e={^)-{e)\ 

Figs. [T] and [5] show the initial eccentricity in depen- 
dence of the beam energy for different centralities and 
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FIG. 1: (Color online) Initial eccentricity erp (full lines) and 
e P art (dotted lines) at time t s tart for Au+Au/Pb+Pb collisions 
as a function of the beam energy. The results are shown for 
four different centrality selections, b < 3.4 fm, 6 = 5 — 9 fm, 
b — 9 — 11 fm and fixed impact parameter b = 7 fm. The 
standard deviation 8e is depicted as error bars to the mean 
value. 
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FIG. 2: (Color online) Initial eccentricity erp (full lines) and 
e P art (dotted lines) at time tstart for Pb+Pb collisions as a 
function of the impact parameter for _Ei a b = 40^4 GeV (a) 
and -Eiab = 160A GeV (b). The standard deviation Se is 
depicted as error bars to the mean value. 



the impact parameter dependence for two different beam 
energies, respectively. The full lines depict the reaction 
plane eccentricity while the dotted lines refer to the par- 
ticipant eccentricity. In these calculations all particles 
which have suffered at least one interaction are included. 
This is consistent with the particles that are taken into 
account for the hydrodynamic evolution where the spec- 
tators are propagated separately in the cascade. 

The initial eccentricity - mean value and standard de- 
viation (as error bars) - as it is shown in Fig. [1] increases 
with the beam energy for both definitions. For simplicity 
the starting times are held constant with respect to the 
centrality variation. As expected the eccentricity is larger 
for more peripheral collisions. The participant eccentric- 
ity is always larger than the reaction plane eccentricity 
since it is calculated in the frame where the eccentricity 
is the largest possible one. Also the fluctuations of the 
eccentricity grow with increasing beam energy and with 
decreasing centrality of the collision. The fixed impact 
parameter calculation for b = 7 fm leads to very simi- 
lar results as the b = 5 — 9 fm calculation for mid-central 
events. This hints to the fact, that centrality fluctuations 
have only a minor impact on the overall event- by- event 
fluctuations of the initial state. 

Apart from the fact that the initial eccentricity at the 
highest SPS energy (see Fig. [2]) is larger than at = 
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40A GeV the results for both energies are pretty similar. 
The reaction plane eccentricity is zero for small impact 
parameters while the participant eccentricity stays finite 
due to the rotation to the participant plane. Going to 
peripheral collisions the participant eccentricity increases 
almost linearly and the reaction plane eccentricity drops 
down again. The fluctuations grow much more for the 
reaction plane definition beyond an impact parameter of 
b = 10 fm. The most interesting observation from Fig. 
[5] is the small difference between e^p and e part and the 
very moderate fluctuations for mid-central collisions in 
the region of b = 5 — 9 fm. This behaviour is directly 
reflected in the centrality dependent elliptic flow results 
as it is shown in Section Hvl 



III. THE HYBRID APPROACH 

To simulate the dynamic evolution of rclativistic heavy 
ion reactions combined microscopic -(-macroscopic ap- 
proaches have proven to be very successful in the descrip- 
tion of various observables |15l425| . The approach that we 
are using here has recently been developed and is based 
on the UrQMD hadronic transport approach including a 
( 3 + 1 ) - dimensional one fluid ideal hydrodynamic evolu- 
tion [2^, [27J for the hot and dense stage of the reaction 
while the early non-equilibrium stage and the final de- 
coupling is treated in the hadronic cascade [H, [29[ [44[ . 

For the present investigation two in principal different 
setups are employed, but since they are constructed from 
the very same ingredients fair comparisons can be made. 
The first configuration is the integrated approach where 
the whole evolution from the incoming nuclei to final 
state particle distributions is calculated on an event-by- 
event-basis. The second possibility is to stop the UrQMD 
calculation at i s tart and to average over many events at 
this time. These averaged initial conditions are then used 
to calculate the hydrodynamic evolution once and then 
the Cooper-Frye transition and the subsequent hadronic 
rescattering is averaged over many events again to ob- 
tain good statistics for the observables. In this second 
setup the spectators are not taken into account for the 
further evolution. To avoid spoiling the results by this 
difference in the setting we concentrate on observables at 
midrapidity. 

In both setups the particles in the UrQMD initial state 
are mapped to energy, momentum and net baryon den- 
sity distributions via three-dimensional Gaussian distri- 
butions that represent one particle each 30]. Two dif- 
ferent equations of state are used to exemplify the dif- 
ferences due to this external input. One is a hadron gas 
equation of state (HG) with the same degrees of freedom 
as in the UrQMD approach [3l| . The other one is a bag 
model equation of state (BM) including a strong first or- 
der phase transition to the quark gluon plasma with a 
large latent heat [13] ■ To see if fluctuations in the initial 
state affect the result differently for different expansion 
dynamics during the hydrodynamic evolution these two 



extreme cases have been chosen. 

The transition from the hydrodynamic evolution to the 
transport approach when the matter is diluted in the late 
stage is treated as a gradual transition on an approxi- 
mated iso-eigentime hyper-surface (see [32[ for details). 
The final rescatterings and resonance decays are taken 
into account in the hadronic cascade. 



IV. ELLIPTIC FLOW RESULTS 

The elliptic flow, the second coefficient of the Fourier 
expansion of the azimuthal distribution of the particles 
V2 , quantifies the momentum anisotropy in the transverse 
plane f33T - l35j ] . It is a self- quenching effect, since the co- 
ordinate space asymmetry is transformed to a momen- 
tum space anisotropy until the system becomes isotropic. 
The elliptic flow is sensitive to the pressure gradients and 
therefore to the equation of state of the matter [36|, l37| ■ 
It has also been shown that the elliptic flow is very sen- 
sitive to the shear viscosity that is present during the 
evolution. 

In Fig. [3]the centrality dependence for the charged par- 
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FIG. 3: (Color online) Centrality dependence of elliptic flow 
of charged particles at midrapidity (\y\ < 0.5) for Pb+Pb 
collisions at E^b = 40A GeV (a) and £i ab = 160^4 GeV (b). 
The horizontal lines indicate the results for averaged initial 
conditions using the hadron gas EoS (blue full line) and the 
bag model EoS (black dotted line) while the symbols (full 
circles for HG-EoS and open squares for BM-EoS) depict the 
results for the event-by-event calculation. 
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FIG. 4: (Color online) Excitation function of charged par- 
ticle elliptic flow for mid-central Au+Au/Pb+Pb collisions 
in comparison to the experimental data (coloured sym- 
bols) [33 . 140] ]. Results for the hadron gas EoS with av- 
eraged/fluctuating initial conditions are depicted as blue 
full/broken line while the calculations with the bag model EoS 
with averaged/fluctuating initial conditions are represented as 
black dashed/dotted line. 



ticle elliptic flow is shown. The upper plot (a) presents 
results from the hybrid model for Pb+Pb collisions at 
Eub = 40^4 GeV while the lower plot (b) is for the highest 
SPS energy. The horizontal lines depict calculations for 
averaged initial conditions, while the symbols represent 
the full event-by-event-setup. Overall, the differences in 
the integrated elliptic flow between the two different cal- 
culations are small. Within the statistical error bars the 
averaged results are in line with the results including 
fluctuations for the corresponding centrality range. In 
accordance with the results for the impact parameter de- 
pendence of the initial eccentricity the fluctuations get 
larger for peripheral events. The highest elliptic flow val- 
ues are reached for mid-central collisions (b = 5 — 11 fm). 
Therefore, we will refer to the impact parameter range 
from b = 5 — 9 fm for the following considerations. 

Furthermore, calculations for the two different equa- 
tions of state are shown. For the bag model EoS ('BM') 
the elliptic flow is expected to be smaller than for the 
hadron gas EoS ('HG') since the speed of sound during 
the expansion is smaller. This expectation is not ful- 
filled if one looks at the results [38|]. The final elliptic 
flow results for different EoS are very similar and agree 
within the error bars. The softer expansion seems to be 
compensated by the longer lifetime of the system. 

The calculation of the integrated elliptic flow at rnidra- 
pidity in dependence of the beam energy (see Fig. 0]) 
confirms the finding that the results do not depend on 



FIG. 5: (Color online) Excitation function of charged particle 
elliptic flow for mid-central Au+Au/Pb+Pb collisions in com- 
parison to the experimental data (coloured symbols) [39l. l4fj| . 
The initial state elliptic flow is shown as the red dashed line. 
The dotted lines refer to the end of the hydrodynamic evolu- 
tion while the full lines represents full hybrid calculations for 
two different EoS. 



the equation of state or on the initial condition setup. In 
the SPS range the results of the hybrid model calculations 
are in a reasonable agreement with the experimental data 
while at lower energies the generated elliptic flow is too 
high. Here, the inclusion of a mean field in the cascade 
part of the evolution help s to reproduce the data [4l[ . As 
it has been shown in [29j the elliptic flow at SPS energies 
might be affected by the transition from hydrodynam- 
ics to the hadronic transport model and leaves room for 
finite viscosities during the hydrodynamic expansion. 

One might now conclude, that the elliptic flow has to 
be generated in the very early non-equilibrium stage and 
in the final state, if the results are not sensitive to the EoS 
during the hydrodynamic expansion and do not change if 
the overall setup is changed - event- by- event vs. averaged 
conditions. In Fig. [5] the contributions to the final ellip- 
tic flow value from the different stages of the evolution 
are shown. The elliptic flow in the initial state for the 
hydrodynamic expansion at £ s +,art is compatible with zero 
at all beam energies. At early times, there is not enough 
transverse expansion to build up elliptic flow. The dotted 
lines show the elliptic flow directly after the Cooper-Frye 
transition at the end of the hydrodynamic expansion and 
at that time the final elliptic flow is already visible. The 
full lines correspond to the full calculation and indicates 
that only a small amount of elliptic flow is generated in 
the final and dilute stage of the collision. So, one is lead 
to the conclusion that the main part of the elliptic flow is 
indeed generated during the ideal hydrodynamic expan- 
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FIG. 6: (Color online) Transverse momentum dependence of 
elliptic flow of pions at midrapidity in mid-central Pb+Pb col- 
lisions at SPS energies (S lab = 40A GeV, (a) and £i ab = 160A 
GeV, (b)). The full lines depict results for fluctuating initial 
conditions while the dotted lines refer to averaged initial con- 
ditions for two different EoS. 



FIG. 7: (Color online) Transverse momentum dependence of 
elliptic flow of protons at midrapidity in mid-central Pb+Pb 
collisions at SPS energies (£i a b = 40^4 GeV, (a) and i?i a b = 
160A GeV, (b)). The full lines depict results for fluctuat- 
ing initial conditions while the dotted lines refer to averaged 
initial conditions for two different EoS. 



sion. 

Let us now investigate the transverse momentum de- 
pendence of the elliptic flow results. The differential flow 
results might be more sensitive to the different setups. 
Fig. |5] shows the elliptic flow of pions at SPS energies 
for the two different equations of state. For the hadron 
gas equation of state the different setups lead to simi- 
lar results again, while for the bag model the averaged 
initial conditions lead to reduced elliptic flow, especially 
at higher p t . The Bag model equation of state seems to 
be more responding to the initial energy density fluctua- 
tions. If there are fluctuations, the transverse slices that 
freeze-out early and produce high transverse momentum 
particles have a longer time to expand because of the 
higher energy density spots whereas for the averaged con- 
ditions the whole system behaves more smoothly. 

In Fig. [7] the transverse momentum dependence of the 
elliptic flow of protons at SPS energies is shown. The 
protons are interesting because they reflect the dynamics 
of the incoming nucleons and the finite net baryon den- 
sity. For the protons the elliptic flow results are even less 
sensitive to either the EoS or the initial condition averag- 
ing than for the pions. At both energies all calculated V2 
curves are compatible with each other within error bars. 
The comparison to experimental data is postponed un- 



til more reliable results become available (for pions the 
comparison has been published in (38j). 



V. SUMMARY AND CONCLUSIONS 

The effect of initial state fluctuations on the finally 
measurable elliptic flow has been studied within a (3+l)d 
Boltzmann + hydrodynamics transport approach. The 
elliptic flow has been calculated as a function of the im- 
pact parameter, the beam energy and the transverse mo- 
mentum for two different equations of state and for av- 
eraged initial conditions or a full event-by-event setup. 
These investigations allow for the conclusion that in mid- 
central (b — 5 — 9 fm) heavy ion collisions the final elliptic 
flow reaches its maximum and the fluctuations due to the 
initial state fluctuations are at a minimum. It has been 
confirmed that most of the vi is build up during the hy- 
drodynamic stage of the evolution. The final integrated 
and differential elliptic flow for charged particles at SPS 
energies seems to be mostly sensitive to viscosity and not 
so much on the equation of state. Therefore, the use of 
averaged initial profiles docs not contribute to the uncer- 
tainties for the extraction of transport properties of hot 
and dense QCD matter based on viscous hydrodynamic 
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calculations, while other ambiguities e.g. the contribu- 
tion of the bulk viscosity [42| and different treatments of 
relaxation times for a multi-component system [43| re- 
main under debate. Still the question if this statement is 
also valid at higher RHIC energies arises, since our ap- 
proach in its present form can only be applied up to top 
SPS energies. 
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